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ABSTRACT 

During its yearlong outburst in 1975-76, the transient source A0620-00 reached 
an intensity of 50 Crab, an all-time record for any X-ray binary. The source has been 
quiescent since. We recently determined accurate values for the black hole mass, 
orbital inclination angle and distance. Building on these results, we have measured the 
radius of the inner edge of the accretion disk around the black hole primary by fitting 
its thermal continuum spectrum to our version of the relativistic Novikov-Thorne thin- 
disk model. We have thereby estimated the spin of the black hole. Although our spin 
estimate depends on a single high quality spectrum, which was obtained in 1975 by 
OSO-8, we are confident of our result because of the consistent values of the inner-disk 
radius that we have obtained for hundreds of observations of other sources: H1743- 
322, XTE J 1550-564, and notably LMC X-3. We have determined the dimensionless 
spin parameter of the black hole to be a* = 0.12 ±0.19, with a* < 0.49 and a* > -0.59 
at the 3a level of confidence. This result takes into account all sources of observational 
and model-parameter uncertainties. Despite the low spin, the intensity and properties 
of the radio counterpart, both in outburst and quiescence, attest to the presence of a 
strong jet. If jets are driven by black hole spin, then current models indicate that jet 
power should be a steeply increasing function of a*. Consequently, the low spin of 
A0620-00 suggests that its jet may be disk-driven. 

Subject headings: accretion, accretion disks - binaries: individual (A0620-00) - black 
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1. Introduction 



A0620-00 is the prototype soft X-ray transient (SXT), which is an eruptive type of X-ray 
binary. For several days in 1975 the flux at Earth from this source was greater than the combined 
total flux of all the other Galactic X-ray binaries, including Sco X-l. A decade later, a dynamical 
study of its quiescent optica l counterpart (V616 Mon) led to the discovery of the first black hole 



(BH) primary in an SXT (IMcClintoc 



c & Remillard 



1986). A0620-00 is one of eight, similar 



short-period BH SXTs (P ovb < 12 hr; iRemillard & McClintockl 120061) . In the optical band, it is 
the best-studied of these systems bec ause its counterpart is bright (Vquiescent = 18.3) and close: D ■■ 
1.06±0.12kpc dCantrell et alJboiob . 



The radial velocity amplitude of the secon dary star is firmly e stablished and has now been 
determined to the remarkable precision of 0.1% (INeilsen et al.l l2008). However, a reliable estimate 
of the BH mass M, which depends on a robust determination of the orbital inclination angle i, has 
remained elusive. Several measurements of i have bee n made by modelin g the ellipsoidal variabil- 
ity of the secondary, but they have been inconsistent (INeilsen et al.ll2008|) . The determination of 
i is complicated by a variable and phase-dependent component of light from the accretion disk. 
Recently, however, a comprehensive analysis of all of the available light curve data (32 data sets 
spanning 30 years) p oints to consistent v alues of inclination and BH mass: i = 51.°0±0.°9 and 
M = 6.61 ± 0.25 M Q dCantrell et alJboioh . 

Our group has published spin estim ates for five stellar- mass BHs ( Shafee et al.l2006 : McClintock et al 



20061 : lliu et alj|2008c iGou et al.ll2009|) . Meanwhile, the spins of several stellar -mass BHs have 
also been obtained by modeling the profile of the Fe K line (see [Miller et al.ll2009L and references 
therein). The spins we find are all quite high, with values of the spin parameter a* in the range 0.7 
to > 0.98. The dimensionless quantity a* = cJ/GM 2 with \aj < 1, where M and J are respectively 
the BH mass and angular momentum (IShapiro & Teukolskyll 19831) . We use the X-ray continuum- 
fitting method, which was pioneered by IZhang et al.1 ( 



1997). Our spin estimates are based on our 



version of the Novikov-Thorne thin accretion disk model (Li et al 



that contain a dominant thermal component (|Steiner et al 



ment of spectral hardening jDavis et al. 2005 : Davis & Hubenv 20061) . We only consider spectra 



2005) and an advanced treat- 



2009a|) and f or which the Eddington 



scaled bolometric disk luminosity is moderate, / = L bo i(a*,M)/L Edd < 0.3 (M cClintock et al.ll20 06). 



Fo r the continuum- fitting method to succeed it is essential to have accurate values of M, i and 
D (e.g.. lGou et al.ll2009|) . such as those reported above for A0620-00. Herein, we use these input 
data and the only suitable, extant X-ray spectrum of the source in order to estimate the spin of the 
BH primary. Although our spin measurement is based on a single spectrum, we have considerable 
confidence in our result for the following two reasons: (1) In our earlier work, we have found 
that our measurements of spin, or equivalently the dimensionless radius of the innermost stable 
circular orbit r isco = Ri SC o/(GM /c 2 )), for a particular source are remarkably consistent over years 
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or decades (e.g jGou et al.ll2009l : ISteiner et al]|2009all2010|) . Most compelling is our recent study 
of the persistent source LMC X-3: Hundreds of observations obtained over a span of 26 years 
by eight different missions give values of r isco , and hence a*, that are consistent within a few 
percent (Steiner et al. 2010). One important conclusion from this study is that a single high quality 
spectrum is a good proxy for a large collection of spectra. (2) The OSO-8 spectrum in question is 
a very high quality spectrum for the determination of spin via the continuum-fitting method: It is a 
remarkably pure thermal spectrum that is almost completely free of the effects of Comptonization 
(Section 3), and it was obtained using a stable, advanced proportional-counter detector that was 
the forerunner of the RXTE PC A detectors. 



2. Data Selection and Reduction 



A0620-00 was observed using the Godd ard Cosmic X-ray Spectroscopy Experiment (GCXSE) 
aboard OSO-8 (|Serlemitsos et al.ll 1976L 1 1977b . The detector employed, namely the C Detector, is a 
sealed multiwire xenon-methane proportional counter with a net effective area of 237 cm 2 , which 
is fitted with a 5.°1 (FWHM) circular collimator. The observation commenced on 1975 September 
29 at 1 1:58 UT and c ontinued for 3.0 days during which the intensity of the source was ~ 1 1 Crab 
( Matilsky et al.|[l976 ). The average collimator-corrected raw count rate was w 1650 counts s" 1 and 
varied by rs 20% over the course of the observations. The gain and detector resolution (18% at 
6 keV) were deter mined using an on-board 241 Am source, and the energy calibration is good to a 
precision of ~ 1% (jMushotzky et al.lll978|) . 



We first ran the FTOOL osofindfast to determine the "good observing days" and then down- 
loaded the appropriate raw PHA data. These data had already been filtered to eliminate times of 
high background during passage through the South Atlantic Anomaly and times of Earth occulta- 
tion. Attending only to the data from Detector C, we ran the tool osopha to extract the 63-channel 
spectrum (spanning 2-60 keV), and we added a 2% systematic error to all the channels; osopha 
automatically corrected the counting rates for dead time by the factor 0.704. The effective source 
exposure time for the observation is 21.3 ks. We then used the tools osofindfast and osopha to 
extract and examine background spectra from several locations. Our results are insensitive to the 
choice of the background region, and we finally settled on a region centered at I = 201.°9 and 
b = -28.°2 and an effective background observation time of 3.35 ks. We computed a correcti on to 
the response of the detector b y comparing the power -law spectrum of iToor & Sewardl (|1974|) . our 
standard reference spectrum (IMcClintock et al.ll2006|) . to the parameters derived by analyzing a 
spectrum of the Crab that was o btained just 12 day s prior to the obser vation of A062 0-00. Follow- 
ing a new correction procedure (ISteiner et aU2010|) . we computed a lToor & Sewardl normalization 
coefficient / TS = 0.792 and a slope difference Ar TS = -0. 148, where / T s is the ratio of the observed 
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normalization to that of llbor & Seward and Ar TS is the difference between the observed value of 
the photon index and that of the reference spectrum. This correction is applied in all the analysis 
work below via a custom XSPEC multiplicative model. 



3. Analysis and Results 



Preliminary, nonrelativistic analysis: All of the da ta analysis and model fitting throughout 
this p aper were performed using XSPEC version 12.6 (|Arnaudlll996|) . As in our earlier work 
(e.g., iGou et al.1 12009|) . we first make an assessment of the data by analyzing it with the non- 
relativistic disk model DISKB B. For modeling the weak Compton component of emission we use 
our convolution model SIMPL (Steiner et al. I l2009bh . which far outper forms the standard power-law 
mode l (powerlaw) with its troublesome divergence at low energies (|Yao et al. boost Steiner et al. 



2009al) . The parameters of SIMPL and powerlaw are similar; their principal parameter, the pho- 
ton index Y is identical. However, the normalization parameter for SIMPL is the scattering fraction 
- the fraction / sc of the seed photons that are scattered into the power-law tail - rather than the 
photon flux. 

Thus, the model we first em ploy is TBABS(SlM PLCjjDDlSKBB), where TBABS is a widely-used 
model o f low-energy absorption (|Wilms et al.ll2000|) . We fitted the data over the energy range 2.2- 
17 keV (|Mushotzky et al.lll978|) . Because of the detector's limited low-energy response, we are 
unable to fit for the hydrogen column density A%, which we estimate from published values of 
the reddening. Based on the five refereed papers known to us, we find 0.25 < E(B — V) < 0.45 
with a most frequently-cited value of 0.35 (|Wu et al.lll983|) . We adopt th e value E(B - V) = 0.35 ± 
0.05. Assuming A v / E(B-V) = 3. 1 and N H /A W = 2.0 x 10 21 mag" 1 cm" 2 dPredehl & Schmittlll995 : 
Giiver & Ozelll2009|) . we obtain our estimate of the column density: _/V H = 2.2 x 10 21 cm" 2 . 



The fit to the PHA spectrum is good (xl = 1-03), and the temperature is precisely deter- 
mined: kT = 0.700 ± 0.004 keV. Moreover, the parameters of the Comptonized emission are well 
determined: T = 3.81 ± 0.3 7 and fcr = 0.0 08 ± 0.003. The best-fit value of the scattering frac- 
tion is low, 0.8% (compare IGou et al.ll2009l and lSteiner et al.ll2009a|) . and it is even lower for the 
relativistic model (0.6%; see be low). Thus, this spectrum is thermal dominant in the extreme 
(|Remillard & McClintockl |2006|) . Because the Compton component is faint the fitted values of 
both kT and a* depend very weakly on how one models this nonthermal emission. 



We are fortunate that in porting the OSO-8 GCXSE software/data to the HEASARC that 
precisely the same spectrum of A0620-00 we consider was used to illustrate how one reduces and 
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analyzes OSO-8 data_|. The model employed in this example is phabs(DISKBB+POWERLAW), 
and the fixed column density is ~ 35% higher, 7V H = 3.0 x 10 21 cm" 2 . Using this model, we obtain 
precisely the same temperature reported in the example, namely kT = 0.70 keV, which confirms 
the correctness of our reduction/analysis procedures. 

Relativistic analysis: We turn now to the analysis of the data using our fully relativistic accre- 
tion disk model KERRBB2, wh ich includes self-irradiation of the disk ("returning radiation") and 
limb darkening jLi et al. 2005 ). The effects of spectral hardening are incorporated into the basic 



model kerrbb via a pair of look-up tables for the hardening fact or / correspondi ng to two repre 



sentative values of the viscosity parameter: a = 0.01 and 0.1 (see lGou et al 



2009). The entries in 



this table were computed using a second relativistic disk model BHSPEC jDavis et aL 2005 ). We 



refer to the model kerrbb plus this table/subroutine as KERRBB2. The model kerrbb 2 has just 
two fit par ameters, namely the BH spin parameter a* and the mass accretion rate M. For further 



details see lMcClintock et al.l (2006). 



We now analyze the data in exactly the same manner as before, except that we replace diskbb 
by kerrbb2: tbabs(SIMPL®kerrbb2). We fix the column density at N H = 2.2 x 10 21 cm" 2 , and 
we have four free parameters: the dimensionless spin parameter a*, the mass accretion rate M, 
the photon index of the high energy component T, and the scattering fraction / S c- We fitted the 
spectrum with the input parameters fixed at their baseline values (see footnote to Table [T]). The 
normalization was fixed at unity (as appropriate when M, i and D are held fixed). We included the 
effects of limb darkening (lflag =1) and returning radiation (rflag =1), and we set the torque at 
the inner boundary of the disk to zero (rj = 0). The fitted results are presented in the first row of 
TableQ] Compared to the nonrelativistic fit, the fit here is even better (xl = 0.74). The luminosity 
is moderate (I « 0.1) and easily meets our selection criterion I < 0.3 (Section 1). Of prime interest, 
we find a low and precise value of the spin parameter, a* = 0. 135 ± 0.029. 

The unfolded photon spectrum and fit residuals are shown in Figure Q] The best-fit value of 
the scattered fraction is very low, only 0.6% (TableQ]). The second and third rows of Table Q] show 
that the uncertainty in 7V H has a quite modest effect on the value of the spin parameter, shifting the 
best-fit value of a* by < 0.22a. The last two rows of Tabled] further show that fixing the power-law 
index T at ±lcr from its best-fit value of 3.53 also has a modest effect on a*, shifting its best-fit 
value by at most « 0.55a. In fact, relative to the total uncertainty in a* these small shifts are over 
ten times smaller still, i.e., ^O.la, as we now show. 

Comprehensive error analysis: The statistical uncertainty in a* is small, and other sources 
of error dominate. Ignoring for now uncertainties in the theoretical model (see Section 4), we 
consider the effects of uncertainties in (1) the input parameters M, i and D, (2) the column density 



1 http://heasarc.gsfc.nasa.gov/docs/oso8/software/oso8_example3.html 
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Nu, (3) the viscosity parameter a (Davis et al.ll2005l) . and (4) the metallicity Z. We first note that, 
as in the case of LMC X-l JGou et alJboolSk the uncertainty in the metallicity of the disk gas is 



negligible compared to the statistical error : We find that a* changes by only 0.34a when we vary 
the metallicity from our default value of Z = 1 (solar) to Z = 0. 1. 

In order to determine the error in a*, we performed an analysis that considers at once the above 
items (1) and (2), i.e., the combined uncertainties in M, i, D and N K - We fix the viscosity parameter 
to its baseline value, a = 0.01. In order to determine the error in a*, we performed Monte Carlo 
simulations assuming that the uncertainties in the four parameters are normally and independently 
distributed. Specifically, we (1) generated 3000 parameter sets for M, i, D and 7V H ; (2) computed 
for each set the look-up table for the spectral hardening factor / using the model BHSPEC; and (3) 
using these /-tables, obtained a* by fitting our model to the spectrum. The results for the 3000 
simulation runs are shown in Figure [2] The uncertainty in D (panel c) dominates the error in a*, 
with the uncertainty in M having less than half as much effect, while the uncertainties in i and 
Afa are smaller still. The corresponding histogram for the spin displacements in a* about the most 
probable value is shown in Figure[3]by the light solid line. 

The analysis above is for our baseline value of a = 0.01. Following precisely the same pro- 
cedures, we also performed the Monte Carlo analysis for a = 0. 1, and the result is the dashed-line 
histogram shown in Figure[3l The summation of these two histograms results in the large histogram 
(heavy solid line). The combined distribution, which corresponds to a total of 6000 simulations, 
has a median spin value of 0.12 and implies a la error of (-0.19, +0.19). Thus, considering all 
significant observational and model-parameter uncertainties, we arrive at our final result for the 
spin of A0620-00: a* = 0.12 ±0.19 (la). The corresponding radius of the inner disk, which is 
uncertain by 12%, is just 7% less than the Schwarzschild value of 6GM/c 2 . 



4. Discussion 



Any measurement of BH spin is only as good as the theoretical model behind it. The continuum 
fitt ing method assumes that th e radial profile of the disk L(R) is given by the analytical form derived 
by iNovikov & Thornd (| 19731 . NT). Because any serious error in the NT model ( e.g., the assump- 
tion of vanishing torque at the ISCO) will lead to large systematic error s in the derived BH spin 
values, we have mounted a major effort to scrutinize the NT model. In IShafee et al.l ( 2008b . we 
reported a 3D GRMHD simulation of a thin disk (H/R ~ 0.06) around a nonspinning BH. We 
showed that the angular momentum profile matches the NT prediction to within ~ 2%, indicating 
that any magnetic coupling across the ISCO is we ak, and we estimated th at the a dditional disk lu- 
minosity due to this weak coupling is only ~ 4%. iNoble & Krolikl (|2009|) and lNoble et al.1 (120 10|) 
have carried out simulations of similarly thin disks, but they conclude that deviations from the NT 
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model are much larger, w ith the specific angular momentum deviating by up to 15%. In a recent 
paper (IPenna et aU2010|) . we report simulations corresponding to a variety of disk thicknesses and 
BH spins. For thin disks with H/R < 0.07 we agree with Shafee et al.'s conclusion that the NT 
model provides an accurate description of both the angular momentum profile and the luminosity 
(see Figures 14 and 15 in Penna et al.). We suggest that the contrary results obtained by Noble 
et al. are because (i) they used for the initial magnetic field in their simulations a topology with 
long-range radial coherence, which we argue is inappropriate for a thin disk, and (ii) they included 
the contribution of the corona, even though the corona is generally believed not to participate in 
the optically thick thermal emission we model in the continuum-fitting method. 

The inner X-ray- emitting portion of a thi n accretion disk is presumably aligned with the spin 
axis of the BH (e.g., LLodato & Pringldl2007|) . In determining the spin of A0620-00 and most 
sources, we assume that the BH's spin is closely aligned with the orbit vector; a misalignment 
of more than several degrees would significantly affect our results (e.g., see Figured). There is 
prese ntly no good evidence for signi ficant misalignments despite two often-cited cases (see Section 
2.2 in Narayan & McClintockl2005 ). A recent population-synthesis simulation study indicates tha t 
the majority of BH binaries have relatively small misalignment angles (< 10°; lFragos et al. 2010 ). 
The approved NASA GEMS polarimetry mission, which is scheduled for a 2014 launch, is ex- 
pected to soon pr ovide the capability to determin e directly the degree of alignment to an accuracy 
of a few degrees Li et alJboOft llCaaret et alJl2009b . 



Based on modeling the 1975-76 X-ray and optical light curves. ISuleimanov et al.1 (|2008l) find 
a low and consistent value of spin, a* ~ 0.1, for our mass of M = 6.6 M and for a fixed value of 
f =U (see their Figure 7). However, for this model their fitted value of the viscosity parameter, 
a ~ 0.6, is significantly higher than the higher fiducial value of a = 0.1 that we have adopted 
(Section 3; see discussion in Secti on 5.3 in Gou et a l. 2009) or the values of a ~ 0. 1 -0.4 based on 
the "best observational evidence" (IKing et al. l2007h . We are unable to estimate a* for larger values 
of a because the required bhspec table models (Section 3) do not exist. However, it is clear (e.g., 
Figure 3) that larger values of a will decrease our estimate of a*. 



Recent GRMHD simulations by Fragile! (I2009|) show that luminous, geometrically-thick disks 
with / > 0.6 (corresponding to H/R > 0.2; McClintock et al. 2006) can falsely appear to host BHs 
with low spins if the inner disk is significantly tilted with respect to the BH spin axis. It is highly 
unlikely that this effect can explain the low spin of A0620-00 because the observed disk luminosity 
is low, / ~ 0.1 (Table 1), and the disk is correspondingly thin H/R ~ 0.04. In this case, as Fragile 
notes, one does not expect the inferred value of spin to be suppressed by the tilted-disk effect that 
he describes. 

Despite its faint corona and low spin, A0620-00 was a bright transient radio source in out- 
burst, decaying from ~ 200 mJy to < 10 mJy in the period from about 12 to 24 days following its 
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discovery on 3 August 1975. A reanalysis of these data by lKuulkers et al.1 (|l999|) indicates multiple 
jet ejections of initially optically-thick components. The authors find that the source was extended 
on arcsec scales, and they infer a relativistic expansion velocity. Even in quiescence, the radio 
source was dete cted at a level of ~ 0.05 mJy, indicating the presence of a partially self-absorbed 
synchrotron jet (|Gallo et al.ll2006l) . If j ets are powered by BH spin, then jet power is likely to in- 
crease dramatically with increasing a* dMcKinnevll2005l : iTchekhovskoy et al.ll2010|) . The expected 
dependence i s so steep that, for a* < 0.4, the jet receives more power from the accretion disk than 
from the BH (IMcKinneyll2005l) . Therefore, given the low spin of A0620-00, it would appear that 
the jet inferred in this source was probably driven by the accretion disk, not the BH. In closing, 
we note that a statistical study by Fender et al.l (I2010|) . although based on data of uneven quality, 
found no evidence that BH spin powers jets. 
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Table 1. kerrbb2 fit results for A0620-00 i 







r 


fsc 


a* 


M c 


Z d 




1 


2.2 


3.53 ±0.38 


0.0062 ±0.0030 


0.135 ±0.029 


1.62 ±0.06 


0.109 


0.74 


2 


1.9 e 


3.51 ±0.38 


0.0060 ±0.0029 


0.141 ±0.028 


1.60 ±0.06 


0.109 


0.74 


3 


2.5 e 


3.55±0.39 


0.0064 ±0.0032 


0.126 ±0.029 


1.63 ±0.06 


0.109 


0.74 


4 


2.2 


3.91 f 


0.0099 ±0.0005 


0.116±0.019 


1.65 ±0.04 


0.110 


0.75 


5 


2.2 


3.15 f 


0.0038 ±0.0002 


0.155 ±0.025 


1.58 ±0.05 


0.108 


0.76 



a M = 6.61 Mq, i = 51.°0, D = 1.06 kpc, and a = 0.01. 

b Column density is in units of 10 21 cm" 2 . 

c Mass accretion rate is in units of 10 18 g sec" 1 . 

d Z = L bo i(a*,M)/L EDD (see Section 1). 

Corresponds to E(B - V) = 0.35 ± 0.05 (see Section 3). 

Corresponds to T = 3.53 ± 0.38 from Row 1. 



j i 

5 10 

Energy(keV) 

Fig. 1 . — Unfolded OSO-8 X-ray spectrum of A0620-00. Top: The histogram shows the model 
fitted to the data with the dashed line representing the thermal component. Bottom: The residuals 
(data minus model in units of xl)- The data for the bin centered at E 12.7 keV were corrupted 
and are ignored (see footnote 1 in the text). 
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Fig. 2. — Effect on the spin parameter a* of varying M, i, D and Af H for the case a = 0.01. (a) 
The upper panel shows a normal distribution for the BH mass M and the lower panel shows a* 
versus mass M for 3000 sets of parameters drawn at random. The central filled circle indicates our 
estimate of the spin a#o=0.135^ obtained from these simulations, (b-d) Same as panel a except 
now for the parameters of inclination angle, distance and column density, respectively. 
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Fig. 3. — Combined error analysis that considers both of our fiducial values for the viscosity 
parameter. The thin solid line is for a = 0.01 and the dashed line is for a = 0.1. The sum of 
these two smaller histograms forms the large histogram. The vertical lines pertain to this latter 
histogram: The vertical solid line indicates the median value of the spin determined by these 
simulations: a* = 0.12; the two dashed lines enclose 68.27% (la) of the spin values centered on 
the solid line and imply an observational uncertainty of (-0. 19, +0. 19); the two dotted lines enclose 
95.45% (2a) and imply an uncertainty of (-0.39, +0.34). The upper and lower limits on the spin 
at the 99.73% (3a) level of confidence are respectively a* < 0.49 and a*> -0.59. 



